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Abstract 

Short-time dynamic scaling behavior of the 3D ±J Ising spin glass is 
studied by Monte Carlo methods. Starting the replicas with independent 
initial configurations with a small pseudo magnetization, the dynamic evo- 
lution of the overlap q(t) between two replicas is measured. The initial 
increase of the overlap q(t) is observed and the corresponding exponent 9' 
is obtained. Prom the scaling relation A = dj z — 9' , the dynamic exponent 
z is estimated. 

PACS: 64.60.Ht, 75.10.Nr, 02.70.Lq, 82.20.Mj 



For systems at criticality, the spatial correlation length and the correlation 
time are divergent. Relevant physical observables present singular behavior, char- 
acterized by a power law and scaling. Traditionally, it was believed that universal 
scaling behavior only emerges in equilibrium or in the long-time regime of the dy- 
namic evolution. However, recent works have revealed that the scaling behavior 
already emerges at the macroscopic early stage of the dynamic evolution, after a 
microscopic time scale t mic [p]]. One typical example is that a magnetic system 
initially at high temperature with a small magnetization is suddenly quenched to 
the critical temperature (without any external fields) and then released to the dy- 
namic evolution of model A . A new critical exponent xq has been introduced 
to describe the scaling dimension of the initial magnetization. It was found that 
at the early stage of the dynamic evolution, the magnetization obeys a power law 

M(t)~m t e ', (1) 

where is the initial magnetization and the dynamic exponent 9' is related to 
Xq by 9' = (x — (3/v)/ z Q. For a large variety of systems, the exponent 9' is found 
to be positive, i.e. the magnetization undergoes an initial increase |], f| |§, ||. 
This fact makes the effect of the initial magnetization m very prominent. 

Another observable which exhibits the critical exponent 9 is the autocorrela- 
tion A(t) =< L~ d J2i Si(0)Si(t) >. Here Si(t) is a spin variable. At the critical 
point, A(t) evolves with a power law || 

A(t) ~ t~\ A = - - 6'. (2) 

Here the initial magnetization is set to uiq = 0. The exponent A of the autocor- 
relation still relates to Xq since the exponent Xq is not only the scaling dimension 
of the global initial magnetization, but also of the local magnetization density 

ill- 

In the past years, the dynamic behavior of spin glasses is one of the impor- 
tant subjects in statistical physics. The aging phenomena, i.e. the dependence 
of dynamic observables on the temperature and field history of the samples, im- 
plies that nonequilibrium nature of the experimental situation is unavoidable in 
spin glasses JJJ. Therefore, the short-time dynamics of spin glasses is of special 
importance. The dynamical exponent A of the autocorrelation A(t) has been nu- 
merically measured by several authors [§, |9], |l(| . Up to now, however, the critical 
initial increase of the order parameter has not been investigated and correspond- 
ingly the exponent 9' has not been obtained directly. The difficulty here comes 
from the fact that the standard order parameter in spin glasses is rather different 
from that of simple spin systems with a second order phase transition. 

recent literature on short-time critical dynamics, the exponent 9' is now often denoted 
as 9. To avoid the confusion with the stiffness exponent 9 in spin glasses, however, we recover 
the notation 9' as used in Ref. IpJ. 
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In this letter we report our novel approach to the short-time dynamics of the 
three-dimensional ±J Ising spin glass. For the first time, we measure the new 
exponent 9'. Using the scaling relation A = dj z — 9' and taking the A from the 
recent reference [[| as input, we estimate the dynamic exponent z. 

The system under consideration is the 3D Edwards- Anderson model described 
by the Hamiltonian 

H = - JijSiSj (3) 

<ij> 

where the spins Si = ±1 are located on the sites of a L 3 cubic lattice with periodic 
boundary condition and the sum is over the nearest neighbors. The interaction 
Jij take randomly +1 or —1 with probability 1/2. 

For spin glasses, the measurement of the new exponent 9' is not straightfor- 
ward. Here the magnetization is not the order parameter and does not present 
the behavior of Eq. ([!]). The standard order parameter of spin glasses is a kind 
of quenched average of the square of the local magnetization ||TT| . It is not clear 
what would be the corresponding behavior of Eq. ([J) . 

Recently, a pseudo-magnetization has been introduced, which may be a candi- 
date of an alternative order parameter of spin glasses. The pseudo-magnetization 
is rather similar to the standard magnetization in simple spin systems, e.g. the 
Ising model. The general idea is that if we can find out a ground state of the sys- 
tem (for fixed couplings), then we define the projection of the spin configurations 
onto the ground state as the pseudo-magnetization 

= ^E<$(*)s?, (4) 

where {Sf} denotes the ground state configuration. In case of the Ising model, a 
ground state is a configuration with all spins in a same direction, i.e., the pseudo 
magnetization is just the standard magnetization. The pseudo magnetization has 
recently been applied to the two-dimensional spin glass to determine the critical 
exponents. The results are promising [|12| . Therefore, we expect that starting 



from a random initial state with small pseudo-magnetization mo, the dynamic 
evolution of this seudo magnetization will obey the power law in Eq. ([!]). 

In this paper, we are interested in three dimensional spin glasses. However, 
for a three-dimensional spin glass, the search for the ground state is NP-Hard, 
i.e. there is no polynomial algorithm for it. For systems with relatively large 
lattice size, it is not possible to find out the exact ground states. Instead, we can 
only obtain approximate ground states. This approximation could cause severe 
correction to the dynamic scaling behavior of m(t). To reduce this effect we intro- 
duce the technique of replicas. For a sample distribution {Jij}, we first find out 
an approximate ground state, then define two replicas by starting from two inde- 
pendent initial configurations, but both with a small initial pseudo-magnetization 
mo- We update these two replicas with the Metropolis algorithm, and measure 
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the dynamic evolution of the overlap 



<K0 = ;pE#(f)S?(t). (5) 

i 

Here Sj and Sf denote the spins for the two replicas. If the projection of the spin 
configuration of each replica on the ground state increases with a power law in 
Eq. (H), it is a natural assumption that the mutual projection q{t) of these two 
configurations should increase with a power law 

q(t)~ml t 26 '. (6) 

In order to verify this assumption, we have performed simulations for the 
two-dimensional Ising model. Starting from two replicas with independent initial 
configurations, both with a small initial magnetization m = 0.02 on the lattice 
L = 128, we update these replicas at the critical temperature with Metropolis 
algorithm until 200 Monte Carlo steps, and measure the dynamic evolution of 
the overlap q(t). Average was taken over 200000 couples of independent initial 
configurations. In Fig. [I], the time evolution of q{t) is plotted on double-log 
scale. After about 10 Monte Carlo steps, the curve shows a good power law 
behavior. From the slope in the time interval [10, 200] we estimate 26' = 0.391(4). 
Here the error is estimated by dividing the samples into three subsamples. This 
result agrees well with the direct measurement of 9' for the Metropolis algorithm 
reported in Ref. IES, 9' = 0.197(1). 



To proceed the simulations for the spin glass we need to search for an ap- 
proximate ground state for each sample distribution {J^}. We employ a simple 



and efficient algorithm proposed in Ref. |T^]. On a lattice L = 16, it takes 
about 1000 second on a 400 MHz ALPHA station to obtain an average energy 
per site e = —1.7829(3), which should be good enough by referring to the esti- 
mated ground state energy e ro = —1.785708(75) for the 3D ±J spin glass [ |i~5H . 
With the ground states {Sf} at hand, we generate an initial configuration by 
randomly setting the Si to be +Sf or — Sf with the probability (1 + m )/2 or 
(1 — mo)/2, respectively. The configuration generated in this way has an averaged 
initial pseudo magnetization mo- 

Rigorously speaking, the critical exponent 9' is defined in the limit tuq — > 0. 
This requires the initial pseudo magnetization m to be small enough. We have 
performed the simulations on a lattice L = 16 with m = 0.06, 0.04 and 0.02, 
respectively. The results show that there is already no 'finite mo effect' for mo = 
0.02. Besides, by analyzing rough results on lattices L = 8 and L = 12 we find 
that on the lattice L = 16 the finite size effect can be ignored. Therefore, in our 
simulations we mainly set m = 0.02 and L = 16. 

In recent references, the critical temperature for the 3D ± J Ising spin glass is 



usually cited as T c = 1.175(25) [|T^, [T7||. We perform our simulations with three 



different temperatures around this T c , i.e at T = 1.05, T = 1.175 and T = 1.3. 
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For each sample of {Jy} we take an average over 10000 initial configurations, 
and the final result is averaged over 1000 samples of the {.%}. In each run, we 
update the replicas until 500 Monte Carlo steps. 

In Fig. H, the time evolution of q(t) for the three temperatures is plotted 
on log-log scale. At T c = 1.175, q(t) shows a power law increase after about 
30 Monte Carlo steps. From the slope in the time interval [30,500], we obtain 
the exponent 29' = 0.17(1), i.e. the exponent 9' = 0.085(5). At T = 1.3, q{t) 
obviously deviates from the power law behavior. This is a typical behavior in 
a second order phase transition for T > T c . At T = 1.05, q(t) seems to show 
still a power law with a similar exponent 29' = 0.18(1). This could mean that 
the critical temperature is actually between 1.05 and 1.175. Indeed, in a recent 
paper I8|| , T c is measured as 1.11(4). In principle, with extensive simulations 
to a longer time, the critical temperature T c can be located from the short-time 
critical behavior JT^, [J. However, this requires simulations on lager lattices and 
is beyond our present power of computing. 

With the critical exponent 9' in hand, we can test the scaling relation A = 
d/z—9'. In Ref . [[|, A has been estimated to be 3.9(1). This value agrees very well 
with the experiment ]20|. Takinging A and 9' as input, we obtain the dynamic 
exponent z = 6.3(2). This value is in agreement with the results z = 6.1(3) in 
01, z = 5.85(30) in [P and z = 6.0(5) in |2 . 



In conclusion, we have numerically simulated the dynamic relaxation process 
of the 3D ±J Ising spin glass starting from random initial states with small 
pseudo- magnetization m . The initial power law increase of the overlap q(t) 
between two replicas is observed and the critical exponent 9' is measured. The 
scaling relation A = d/z — 9' is confirmed. 
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Figure 1: The time evolution of the q(t) for 2D Ising Model at T c is plotted in 
double-log scale. The lattice size is set as L = 128 and the initial magnetization 
mo for each replica is 0.02. 
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Figure 2: The time evolution of the q(t) for 3D ± Ising spin glass is plotted in 
double-log scale. From the top to the bottom, the temperatures are T = 1.05, 
1.175 and 1.3 respectively. 



7 



